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ABSTRACT 
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Context. Extrasolar planets observation and characterization by high contrast imaging instruments is set to be a very important subject 
in observational astronomy. Dedicated instruments are being developed to achieve this goal with very high efficiency. In particular, 
full spectroscopic characterization of low temperature planetary companions is an extremely important milestone. 
Aims. We present a new data analysis method for long slit spectroscopy (LSS) with coronagraphy, which allows characterization 
of planetary companions of low effective temperature. In a speckle-limited regime, this method allows an accurate estimation and 
subtraction of the scattered starlight, to extract a clean spectrum of the planetary companion. 

Methods. We performed intensive LSS simulations with IDL/CAOS to obtain realistic spectra of low (R = 35) and medium (R = 400) 
resolution in the J, H, and K bands. The simulated spectra were used to test our method and estimate its performance in terms of 
contrast reduction and extracted spectra quality. Our simulations are based on a software package dedicated to the development of 
SPHERE, a second generation instrument for the ESO-VLT. 

Results. Our method allows a contrast reduction of 0.5 to 2.0 magnitudes compared to the coronagraphic observations. For M0 
and GO stars located at 10 pc, we show that it would lead to the characterization of companions with T e ff of 600 K and 900 K 
respectively, at angular separations of 1.0". We also show that errors in the wavelength calibration can produce significant errors in 
the characterization, and must therefore be minimized as much as possible. 

Key words, techniques: spectroscopic - techniques: image processing - methods: data analysis - stars: planetary systems 



1. Introduction 

Detection of exoplanets through indirect methods, such as radial 
velocities or transits has become a subject of intensive research 
that has led to the discovery of more than 270 exoplanets. A fu- 
ture generation of instruments will be able to deliver very high 
contrast coronagraphic images of nearby stars to detect possi- 
ble orbiting companions. These direct detections are particularly 
effective at detecting relatively hot companions at large separa- 
tions from their host stars, making them the perfect complement 
to radial velocity and transit methods that are usually able to de- 
tect planets with short orbiting periods. 

The European project SPHERE (Spectro-Polarimetic High 
contrast imager for Exoplanets REsearch) for the VLT second 
generation of instruments aims to detect exoplanets down to a 
contrast of 10~ 6 at angular separations as small as 0.1" with 
dedicated extreme adaptive optics and coronagraphs ( |Dohlen| 
[eTaLl [2006) 1. One of the 3 science channels of SPHERE wiU 
be equipped with a near infrared camera, IRDIS (InfraRed Dual 
Imaging Spectrograph), that will offer several observing modes 
for exoplanet detection and characterization including long slit 
spectroscopy (LSS). The prime objective of IRDIS is exoplanet 
direct detection with simultaneous differential imaging (SDI), a 
technique that has been extensively described in the literature 
( Racine et al.|1999 Marois et al.||2000"| including possible im- 
provements ( jMaroFs et al. 2006J. It was also used successfully to 
detect cool companions with the NACO instrument on the ESO- 
VLT ( Lenzen et al.|2004) >. This technique relies on the fact that 
the planetary spectrum has deep absorption features due to its 
atmosphere composition (e.g. CH4 in the H band), whereas the 
starlight has a relatively smooth spectrum over the observation 
band. Recording two simultaneous images inside and outside of 



the planet absorption band, and subtracting them allows removal 
of most of the starlight, leaving mainly the planetary light in the 
subtracted image residuals. 

Although very powerful at detection, this technique appears 
less suitable for the characterization of exoplanets. Current mod- 
els predict that giant gas planet atmospheres most probably con- 
tain methane and water ( Allard et al. 2001 , 2003; [Burrows et aT) 
|1997[ |2006| ). The reliance on the related molecular features for 
detection is a primary step in characterization, since it identifies 
planetary objects with atmosphere deprived of these compounds. 
However, a full characterization is not possible with a flux dif- 
ference of only 2 narrow and adjacent bands, which is the in- 
formation provided by the SDI technique: some degeneracy in 
planet parameters is probable if only the methane band in H is 
observed. It can also be difficult to exclude signal contamination 
by a background object without more information. This explains 
why follow-up spectroscopic observations over a significantly 
wide wavelength range are required. 

In IRDIS, this observation mode is achieved with long slit 
spectroscopy at small (R = 35) and medium (R = 400) reso- 
lutions and coronagraphy over J, H and K bands. In this con- 
text, we adapted and improved the data reduction technique first 



developed by Sparks & Ford ( 2002 1 and used by IThatte et al 



(2007) with AO-fed Integral Field Spectroscopy (IFS), to extract 
a full companion spectrum from simulated high-contrast images. 
In contrast to its use with an IFS, LSS provides spatial informa- 
tion in only one dimension (along the slit), making it necessary 
to develop a specific data analysis method. 

Our method is very general and can be used on corona- 
graphic images as well as non-coronagraphic ones. Moreover, 
it is fairly simple to implement, and does not require a large 
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Fig. 1. Low resolution {R = 35) simulated spectrum with coron- 
agraphy in the J, H and K bands of an MO star with a relatively 
bright companion at an angular separation of 1.0". Wavelength 
is along the horizontal axis, and angular separation along the 
vertical axis, starting from the center. The obscured central part 
corresponds to the position of the 0.4" coronagraph mask. The 
speckles form oblique lines in the spectrum because their an- 
gular separation scales linearly with wavelength, whereas the 
planet position is fixed. For clarity, the displaying scale is not 
linear. 



number of operations. It is fully described here in the context of 
SPHERE, since it is intended to be used for LSS data reduction 
on this particular instrument. Section [2] describes the concept 
and practical implementation details, Sect. [3] provides a global 
overview of the simulations performed to simulate realistic LSS 
data, and Sect. |4] details the global performance of the method 
and the influence of a few key parameters. 



2. Data analysis of coronagraphic long slit 
spectroscopy 

2.1. Speckles in long slit spectroscopy images 

The main limitation in high-contrast coronagraphic images orig- 
inates in the speckles (Soum mer et al. |[2Q07) > induced by atmo- 
spheric phase residuals and instrumental static and quasi-static 
aberrations not corrected by adaptive optics (AO). Most of the at- 
mospheric residuals average out with time, producing a smooth 
halo on which the static and quasi-static speckles (sometimes 
called super-speckles in the literature) are superimposed. The co- 
herence time of the super-speckles extends from tens of seconds 
to several minutes (Ma cintosh et al.|2005 Hinkley e t al.|2007 i, 
so it is safe to assume that they do not average during a typical 
integration. In the case of a single image at a given wavelength, 
the typical angular size of a speckle is 9 spec kie = A/D, where A 
is the wavelength of observation and D is the telescope diame- 
ter. The planet is assumed to be an unresolved point source and 
its image has a full width at half-maximum (FWHM) that is ap- 
proximately equal to the diffraction limit. Differentiating a faint 
planet point-spread function (PSF) from a speckle is therefore 



impossible with a single image, and without using a temporal 
sequence or other spectral information. 

The position and size of a speckle in the field of view is 
wavelength-dependent: as the wavelength increases from A\ to 
A2, the FWHM of a single speckle and its angular separation 
from the star both increase by a factor /I2M1 . Following the 
wavelength axis at a fixed angular separation from the star, one 
would measure a strong modulation in the intensity as the Airy 
rings and speckles cross this position (see Sparks & Ford 2002 



Fig. 24). A possible faint companion signal would remain unde- 
tected in the modulations due to the wavelength dependence of 
the speckle pattern. However, a fixed physical object (e.g. a plan- 
etary companion) will not change its position with wavelength: 
only its FWHM will be multiplied by the aforementioned factor. 
This property enables a good determination and subtraction of 
wavelength dependent features at a given position, and allows 
the detection of physical objects in the vicinity of the star. 

In the case of long slit spectroscopy, an infinity of images 
of the slit at increasing wavelengths are dispersed and superim- 
posed to create a spectrum such as that shown in Fig. [T] The 
wavelength-dependent artifacts such as speckles create oblique 
line^jin the spectrum. A bright speckle will see its size increased 
with wavelength, creating a bright line in the spectrum with in- 
creasing width. A physical object creates a straight line at a fixed 
radial separation from the star, and its FWHM increases linearly 
with wavelength. Our method uses these geometrical properties 
to separate the planet and the star spectra, as proposed by Thatte 
et al.| ( 2007) with spectral deconvolution (SD) in the case of IFS 
observations. 



2.2. Method description 
2.2.1. General concept 

In the image representing the spectrum, each column of pixels 
corresponds to a spectral interval AA of a few nanometers, cen- 
tered on a particular wavelength. From now on, we assume that 
a given column of pixel i in the spectrum corresponds to a sin- 
gle wavelength equal to the central wavelength Aj of the spectral 
interval A/l,. This approximation is valid as long as AAj/A/ <K 1, 
which is the case for the wavelength range and spectral interval 
considered in our work (see Sect. [3]) where AAj/A, < 0.015 in 
any case. 

The first step corrects the spectral dependency of the speck- 
les by spatially rescaling each column of pixels by a factor 
or,- = Aq/Ai, where Aq is the shortest spectrum wavelength. Since 
or, is not an integer for most columns, each column was resam- 
pled in the spatial direction to rescale its length from one inte- 
ger number of pixels to another integer number of pixels. The 
number of resolution elements is fixed by the instrumental setup 
(i.e. the platescale) in the spatial direction, and by the dispersive 
element in the spectral direction. In fact, the spectrum is resam- 
pled within a grid ~15 times smaller than the original pixel grid 
before rescaling the columns by the factor a,-. When calculat- 
ing the new integer number of pixels for column i, this spatial 
oversampling allows us to reduce the round-off error introduced 
by the non-integer factor a,. With the approximation mentioned 
in the previous paragraph, neglecting the spectral interval AA, 



1 In the case of a spectrum with constant spectral interval (case con- 
sidered in this work), the wavelength dependent artifacts create straight 
oblique lines. In a more general case where the spectral interval may 
vary (constant resolution for instance), the artifacts follow a curve as a 
function of wavelength. 
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Fig. 2. Same spectrum as Fig. [7] where each column has been 
rescaled according to its corresponding central wavelength to 
correct for the spectral dependence of the speckles. The speckles 
are now along straight horizontal lines, and the planet follows a 
1 1 A law. The displaying scale is the same as in Fig. [T] 



translates into an error of less than 1.5% in a,-. The final rescaled 
spectrum is shown in Fig. [2] 

The rescaling corrects the spectral dependence of each 
speckle, ensuring that they follow horizontal straight lines, while 
the companion spectrum follows a 1 I A law. Speckles are induced 
by random phase errors in the incoming wavefront that scatters 
the starlight. The fraction of light scattered in each speckle com- 
ing from the companion is negligible compared to the scattered 
light of the star, so it is safe to assume that the spectral decompo- 
sition of a speckle reflects only the star spectrum, multiplied by 
atmospheric and instrumental transmission. The only difference 
between one angular position and another is the amplitude of the 
modulation created by the succession of bright and dark speckles 
superimposed on the star halo and the Airy pattern. As a result, 
the star spectrum can be recreated precisely by averaging all the 
lines for which a full spectrum is available. If the spectrum ex- 
tends to an angular separation +p rad from the star, the lines 
containing a complete spectrum are those with a separation of 
+pD/Ai cycles per pupil diameter, where A[ is the longest wave- 
length in the spectrum. In the case of coronagraphic LSS, we 
have to exclude the central lines corresponding to the position of 
the Lyot coronagraphic mask. We finally obtain a linear average 
spectrum of the star, which we refer to as the model spectrum. 

The next step is to estimate the exact amount of starlight at 
each angular separation from the star. To do so, we extract a lin- 
ear spectrum at each angular position at which we fit our model 
spectrum in amplitude using a least-squares estimation. At the 
beginning, both the linear and the model spectra are normalized 
to unity, and we allow our model normalization to vary between 
a factor of 0.5 and 3 times the initial normalization, in steps of 
1%, during the fitting process. This allows us to optimize the fit 
of the model to the linear spectrum. Our least-squares statistic is 
defined to be 




s 2 = 



zf =1 ^ model -/; linear ] 2 

N 



(1) 



Fig. 3. Final spectrum after removing scattered starlight and 
scaling its dimensions back to their original size. Most of the 
starlight has been removed, and the companion spectrum can be 
easily extracted. Some residuals are visible, especially close to 
the center where the star lies. The displaying scale is the same 
as in Fig. [T] 

where f model and £ linear are the model flux and the flux of the 
linear extracted spectrum, at the considered angular separation 
in each wavelength bin z; /3 is the varying normalization applied 
to the model, ranging from 0.5 to 3 in steps of 0.01; and N is 
the total number of wavelength bins. The parameter s 2 therefore 
evaluates the average deviation between the model and the linear 
spectrum. The optimal fit to the data is obtained by finding the 
normalization f3„, that minimizes s 2 . When the appropriate nor- 
malization p„, has been determined, the model is multiplied by 
the product of this factor and the initial factor used to normalize 
the linear spectrum to unity. 

By doing so for every angular separation in the spectrum, 
we construct a "clean" spectrum i.e. a spectrum devoid of light 
from the companion, that we will refer to as the reconstructed 
spectrum from now on. The reconstructed spectrum is subtracted 
from the rescaled spectrum, to remove most of the scattered 
starlight and reveal the spectrum of the companion. Finally, each 
of the resulting spectrum columns needs to be scaled back to its 
original size to be able to recover the companion spectrum along 
a straight horizontal line. It is then possible to extract a final 
spectrum of the companion unaffected by any scattered starlight. 
The final spectrum is illustrated on Fig. [3] 

2.2.2. Discussion and implementation details 

A few important details need to be emphasized to achieve a ro- 
bust implementation of this data analysis method. First of all, 
this method is affected by the same conditions about the bifur- 
cation point defined by Tha tte et al.| p007). The issue is not to 
estimate the starlight uncontaminated by light from the compan- 
ion, which is achieved by averaging the spectrum, but to have 
a sufficient number of points to complete an accurate amplitude 
fitting. To estimate reliably the amount of scattered starlight at a 
given angular position, it is necessary that at least one data ele- 
ment that is uncontaminated by the light of the companion. The 
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more data points used to fit the linear average spectrum, the more 
accurate the estimation of the stellar contribution, and the more 
reliable the characterization of a companion will be. 

Secondly, for very bright companions, a mask is applied to 
the data being fitted to ensure that the light of the companion 
does not contaminate the amplitude fit. If a companion is very 
bright compared to the scattered starlight that we wish to evalu- 
ate, the fitting will be completely biased: the reconstructed spec- 
trum will be overestimated and the companion spectrum contin- 
uum will be underestimated. When subtracting the reconstructed 
spectrum from the observed data, part of the companion light is 
then removed. This issue is important only for very bright com- 
panions, whose position can be easily ascertained from the data 
to create a mask at the appropriate position. The exact size of 
this mask can be deduced from the instrumental PSF. We choose 
to consider a size of 2.5A/D to completely mask the PSF core 



Table 1. Simulated long slit spectroscopy modes. 



(see Sect. 4.2.1 for a more complete analysis). 

Finally, depending on the slit size and spectral range, it is 
important to consider that speckles located at the edge of the slit 
move out of the slit as the wavelength increases. For the slit sizes 
(0.09" and 0.12", see Sect. [3j and wavelength range (J, H and 
K bands) considered, a significantly large number of speckles 
move out of the slit (~50%). The impact of this effect is iden- 
tical at all positions in the slit, so the final influence is negligi- 
ble on data analysis. However it means that the measured model 
spectrum is not equivalent to the true star spectrum. And that dis- 
crepancy between the two increases with wavelength as speckles 
move out of the slit. This implies that precise photometry of the 
star cannot be performed on the model spectrum. 



3. LSS simulations 

Intensive simulations were completed to simulate realistic spec- 
tra for testing the method. We used the complete simulation 
model of the VLT-SPHERE instrument developed for end-to- 
end simulations. This model is a diffractive code written in 
IDL (Interactive Data Language) that is based on the CAOS 
(Code for Adaptive Optics Systems) problem solving environ- 
ment (Carbille t et al.||2004|l with a package developed for the 
SPHERE project ( |Carbillet et al.|2058) . We do not to describe 
in detail the content of this simulation code, but present a global 
overview of the simulated. 

The SPHERE package for CAOS is a diffractive end-to-end 
simulation code that takes into account multiple sources of aber- 
rations, such as atmospheric residuals, AO correction and optical 
aberrations. It consists of separate modules, which simulate dif- 
ferent parts of the instrument i.e. the extreme AO system SAXO, 
the optical common path, and the science modules. A specific 
part of the code is dedicated to LSS simulations with coronagra- 
phy. This code does not simulate temporal variation in the aber- 
rations: only the AO-filtered atmospheric residuals are changed 
during the simulation. The final output of the code are normal- 
ized images of the PSF, the coronagraphic PSF, and slit images 
at different wavelengths. 

A second code was developed to scale these normalized im- 
ages to the correct photometric values according to the instru- 
ment design, the star being observed, the distance, the angular 
separation, and the T eJ y of the companion. In particular, the code 
accounts for the global throughput of the instrument, as well 
as the atmospheric transmission. The code considers OH lines 
but not their variability. The spectra used to model the compan- 
ions in the simulations are the latest synthetic spectra generated 
by the ENS Lyon group with the PHOENIX code ( |Allard et alT] 



Mode 


Band 


Resolution" 


AA 


Slit width 








(nm/pixel) 


(") 


LRS 


JHK 


35 


13.6 


0.12 


MRS 


J 


400 


1.19 


0.09 



Resolution is given at A = 950 nm. 



2001 ). More precisely, we used the models designated DUSTY- 
2000, COND-2002, and SETTL ( jAllard et al.||200T] [2003] and 
private communication). Finally, a realistic amount of noise was 
added depending on the characteristics of the particular next- 
generation AO imaging instruments i.e. detector noise, sky back- 
ground and instrumental thermal background. The final output 
was a two-dimensional spectrum of the star with one or more 
companions at different angular separations. 

We simulated 0.12" and 0.09" wide slits with a 0.20" radius 
Lyot coronagraph at the center, for an 8 meter diameter tele- 
scope. The number of wavelengths and the spectral interval was 
carefully chosen to produce spectra at low and medium resolu- 
tion, refered to as LRS and MRS respectively (see Table[T]i. The 
amount of optical aberrations was set to be ~50 nm RMS be- 
fore the coronagraph, and ~40 nm RMS after the coronagraph. 
The atmosphere was simulated by a set of decorrelated phase 
screens, each corresponding to approximately 80 nm RMS of 
wavefront aberrations after AO correction. We limited our sim- 
ulations to 100 phase screens to produce a smooth star halo. In 
theory, this number is too small to produce a "true" long ex- 
posure. If we assume that atmospheric residuals have a correla- 
tion time of the order of 1 ms and that a typical exposure time 
for a DIT is 10 seconds, then 10,000 decorrelated phase screens 
would be required to produce a long exposure. However, in prac- 
tice using fewer phase screens has proven to give good results 
in various conditions, while significantly reducing the comput- 
ing time of the simulations. The total amount of aberrations is 
~ 1 10 nm RMS, which corresponds to a Strehl ratio of 83% at 
1.6 yum. Although high compared to existing instruments with 
conventional AO (typical Strehl ratio for VLT-NACO in stan- 
dard conditions is between 40% and 60% in K band, see |Clene"t| 
et al.|2004|, these overall performances are realistic for new in- 



struments with high order AO systems and exceptionally high 
optical quality such as VLT-SPHERE. To generate a significant 
amount of data we created 6 independent sets of data in LRS 
with our simulation code by changing the random number gen- 
erator seeds for every set. Each set produced a spectrum with a 
unique realization of the random atmospheric and instrumental 
aberrations. For MRS, only 2 different data sets were generated 
in J band because of the considerable number of wavelengths 
that had to be simulated. 

Photometry for each data set was calculated to be representa- 
tive of one hour exposure time for M0 and GO stars at a distance 
of 10 pc, and companions with effective temperatures ranging 
from 400 K to 2700 K at angular separations of 0.5", 1.0", and 
1.5" (transposing to 5, 10 and 15 A.U. at 10 pc respectively). 
Table |2]lists the atmosphere models used, their effective temper- 
atures, and contrast in the JHK bands around M0 and GO stars. 
Different noise sources were added as residual variance simu- 
late the detection process realistically i.e. photon noise, flat field 
noise (0.1%), and read-out noise (15 e~). Thermal background 
coming both the sky and the instrument were also added: the in- 
strumental thermal background was simulated to be photon noise 
induced by a constant emission with a fixed number in K band: 
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Fig. 4. Contrast curves for a M0 star at 10 pc in J band (left, A — 1.27 yum), H band (middle, A — 1.58 /im), and K band (right, 
A — 2.10 jt/m), for an exposure time of 1 hour at LRS. The hatched area is the part covered by the coronagraph mask. The plot shows 
the non-coronagraphic PSF (plain line), the coronagraphic PSF (dashed line), and the contrast level after applying our method 
(dotted line). The simulated companions with different T e ff are also plotted at 3 angular separations, 0.5", 1.0" and 1.5", with their 
effective temperature and the atmosphere model. 



344 photon sec -1 pixel -1 ; the sky background has been simu- 
lated as a two-dimensional spectrum calculated to match the sky 
emission at ESO-Paranal observatory in Chile (86 x 10~ 16 pho 



Angular separation in arcsec 



ton sec 1 m 2 as 2 at A = 2.2 //m, values calculated from Patat 



2008). The final output of the simulation represents science cal- 
ibrated images that were reduced following the usual steps of 
data reduction. 



4. Results 

4.1. Data analysis performance 

The global performance takes into account two aspects: (a) the 
achievable contrast, and (b) the fidelity of the extracted compan- 
ion spectrum with respect to the original spectrum. We discuss 
our results in terms of these complementary aspects in the fol- 
lowing sections. 

4.1.1. Contrast reduction 

The goal of the data analysis in LSS was to remove the scattered 
starlight and to improve the contrast. If we consider starlight to 
be the main noise source, this operation will improve the signal- 
to-noise ratio (SNR). Figures |4]and|5]present the 5cr detectivity 
at different wavelengths for a M0 star at 10 pc and companions 
of increasing T e ff at 3 angular separations, at LRS and MRS re- 
spectively. The contrast curve was obtained by measuring the 
standard deviation of the corrected spectrum in a window about 
the considered wavelength of a width equal to that of the slit 
image on the detector. All curves were normalized to the non- 
coronagraphic PSF. 

At LRS, we see a clear improvement in the contrast in J and 
H bands between the coronagraphic PSF before and after correc- 
tion. The gain is of the order of 2.5 magnitudes close to the star, 
and 1 magnitude farther out. Companions as cold as 600 K could 
be characterized adequately in these bands very close from the 
star, while for a temperature of 500 K they are detectable at sepa- 
rations larger than 1.0". In K band, where the contrast is already 
favorable, the method is almost inefficient. For the GO star (not 
plotted here), the results are similar in terms of magnitude gain. 
The results are of higher quality in K for the GO star, where a 
gain of 1 magnitude is achieved at separations larger than 0.7". 
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Fig. 5. Contrast curves for a M0 star at 10 pc in J band (A = 
1.27 /zm), for an exposure time of 1 hour at MRS. The hatched 
area is the part covered by the coronagraph mask. The plot shows 
the non-coronagraphic PSF (plain line), the coronagraphic PSF 
(dashed line) and the contrast level after applying our method 
(dotted line). The simulated companions with different T e ff are 
also plotted at 3 angular separations, 0.5", 1.0", and 1.5", with 
their effective temperature and the atmosphere model. 



At MRS the gain is comparable to the LRS case in J band 
with a 2.5 magnitude increase in the contrast close to the star, 
but a constant level is reached at an angular separation of ~0.5". 
For larger separations, the final contrast level remains equal to 
10 . Further simulations show that this is the noise level of the 
sky that sets this limit. Using a thinner slit would decrease the 
sky flux, and then its noise level, probably allowing us to reach 
performances comparable to the LRS at the same angular sepa- 
rations. 

Although calculated for a given instrumental configuration, 
these contrast plots are representative of our method perfor- 
mance because for a given configuration they show a clear im- 
provement in the contrast. The limiting factor here is probably 
the correlated speckle residuals that appear on Fig. [3] These 
residuals come from an inaccurate estimation of the scattered 
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MO a 10 pc + COND 500K 1.0 




2.0 2.5 3.0 




pc + COND 1000K 9 0.5" 




Fig. 6. Example of quality factor for 3 different companions around an M0 star at 10 pc in LRS. In the upper-right part of each plot, 
we indicate the correlation factor p, the discrepancy factor d, and the final quality factor q. The plot shows the spectrum that was 
introduced as input in the data analysis (planet only, plain line), the spectrum before data analysis (star and planet, dotted line), and 
the spectrum extracted as a result of data analysis (dashed line). 



Table 2. Planetary models included in the simulations. 



Model 


Teff 
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Contrast" M0 
(Mag) 
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SETTL 


2400 
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2500 
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Contrast is calculted over the J, H and K bands. 



light during the amplitude fitting of the model spectrum. Since 
the model spectrum is an average, it is likely that small varia- 
tions in the speckle spectrum at each angular separation cannot 
be measured exactly. 



4.1.2. Extracted spectrum quality 

The quality of the data analysis must be quantified to evaluate the 
confidence that can be attributed to a given result. Two different 
aspects of extraction quality need to be monitored: (a) the corre- 
lation and (b) the discrepancy between input and output spectra. 
The former is measured by a simple correlation coefficient 



Ps„s„ = 



E[SiS„]-E[Si]E[S ] 



(2) 



where 5, is the input spectrum, S a is the output spectrum, and 
E [S ] and cr s denote the mean and standard deviation of a spec- 
trum S respectively. The latter aspect is measured by the factor 



d : 



S;,S„ 



\Si-S \ 



Max (Sj) 



(3) 



This factor represents the mean discrepancy between the input 
spectrum and the spectrum that is extracted during analysis. It 
is normalized such that it provides a value between and 1 ; the 
closer to the value is, the higher is the quality of the data anal- 
ysis. The final quality factor is obtained by multiplying the two 




Fig. 7. Example of quality factor in J band for an M0 star at 10 pc 
and a 1000 K companion at 1.0" in MRS. At this resolution, 
spectral features are clearly visible. In the upper-right part of the 
plot is given the correlation factor p, the discrepancy factor d and 
the final quality factor q. The plot shows the spectrum that was 
introduced as input in the data analysis (planet only, plain line), 
the spectrum before data analysis (star and planet, dotted line) 
and the spectrum extracted as a result of data analysis (dashed 
line). 



previous factors to measure both the correlation quality and the 
discrepancy between the input and the output: 



qs„s„ = Ps„s„ (1 -d s „sj- 



(4) 



A value of q = 1.0 denotes a perfect data analysis (i.e. 5, 
and S are identical), while smaller values denote lower levels 
of quality. It is necessary to define a ground value above which 
the data analysis can be considered successful or adequate. 

Figure|6]presents three cases with very different values of the 
quality factor at LRS. The lowest quality case corresponds to a 
500 K companion at 1.0" from an M0 star at 10 pc, for which we 
have a 49.47% quality factor. The most prominent features in the 
input spectrum appear in the extracted spectrum, giving an inter- 
mediate correlation factor of 63.19%, but the mean discrepancy 
between the input and output spectra is quite high, reaching more 
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Table 3. Error on T e jj estimation as a function of quality factor. 



Quality factor 


T e jj error 


(%) 


(K) 


0-9 


1400 


10-19 


800 


20-29 


600 


30-39 


1100 


40-49 


200 


50-59 


300 


60-69 


200 


70-79 


100 


80-89 


100 


90-100 


100 



than 100% in some areas of the spectrum. Photometry analysis 
of such a spectrum would produce spurious results. The second 
case is a 700 K companion at 1 .0" from the same star. The corre- 
lation is high, and the mean discrepancy is below 10%, leading 
to a high quality factor of 80.85%. All the main features of the 
spectrum are reproduced with good accuracy, except for some 
local residuals that have not been completely eliminated during 
data analysis. Finally the most successful case is a 1000 K com- 
panion at 1 .0" from the same star, which has a quality factor of 
96.54%. All features in the spectrum are reproduced with high 
precision, allowing an optimal characterization of the compan- 
ion from its spectrum. The quality factor for a given star clearly 
depends on the companion T e jf and angular separation, i.e. on 
the overall contrast. 

Figure [7] shows a good case of data analysis at MRS around 
a GO star at 10 pc. The discrepancy is low (less than 2%) and the 
correlation almost 100%, providing a final quality factor of more 
than 98%. The dotted line that represents the spectrum before 
data analysis clearly indicates the gain provided by our method: 
all scattered starlight has been removed, and the spectrum shows 
distinct spectral features that were lost in the starlight. 

One of the primary goals of using LSS on planetary compan- 
ions is to determine accurately the T c y/ of the target. We studied 
the error in the T t yy estimation as a function of the quality factor 
value. To estimate this error, we recalculated photometry with 
the COND model from the ENS Lyon group for effective tem- 
peratures ranging from 400 K to 1500 K in steps of 100 K. Each 
of these COND models were introduced into our data analysis 
code for M0 and GO stars at 10 pc and angular separations of 
0.5", 1.0", and 1.5". The output spectra of the data analysis was 
then compared with the complete COND librar^to estimate the 
T e ff from the extracted spectra. Table [3] provides the absolute 
value of the error in T f yy as a function of the quality factor (in 
bins of 10%). For small quality factors (below 40%), the error 
in T e ff is significant. Results with quality factors in the range 
- 40% should therefore not be trusted. For quality factors in the 
range 40 - 70%, the error is ~200 K, which represents an error 
of more than 25% for a companion with T e yy = 700 K. For a ac- 
curate estimation of T e yy , quality factors of more than 70% are 
required. Above this limit, the error is 100 K or lower. To achieve 
yet higher accuracy, a finer model grid in temperature should be 
used. In the following, we choose to set the lower limit for the 
quality factor to be 80%, i.e. we consider that results of quality 
factors smaller than 80% do not provide a sufficiently precise 
estimation of T e ff . 




400 600 800 1000 

Effective Tefferature (K) 



2 The complete library contains models with T e yy from 100 K 
to 3000 K in steps of 100 K. It is available at ftp://ftp.ens- 
lyon.fr/pub/users/CRAL/fallard/ 



Fig. 8. Quality factor for companions around an M0 star at 10 pc 
in LRS with different T f yy at 3 angular separations, 0.5" (plain 
line), 1.0" (dot line) and 1.5" (dashed line). The horizontal 
dashed line marks the 80% level above which the estimation of 
the T e ff can be accurately determined. Error bars represent the 
standard deviation of the quality factor over the 6 independent 
data sets that we have simulated in LRS 



Figure [8] indicates the quality factor values for compan- 
ions with various effective temperatures at 3 angular separations 
around an M0 star at 10 pc in LRS. The plot gives an overview 
of the evolution of the quality factor with effective temperature 
and angular separation. The evolution with T e yy is globally the 
same at each separation: the quality factor slope is very high 
over a small range of effective temperature, going from close 
to to approximately 80%. Above this limit, the slope slowly 
decreases to reach an almost flat regime above 95%. The dif- 
ference with angular separation is the temperature at which the 
80% level is reached: there is a large improvement in contrast 
between 0.5" and 1.5". At 0.5", the 80% level is reached for a 
T e ff of 700 K (contrast of 11.0 magnitudes), while at 1.5" the 
same level is reached for a T e ff of 500 K (contrast of 13.6 mag- 
nitudes). Results are of even higher contrast for the GO star (not 
plotted here), generating wider error bars, in particular for low 
T e ff for which the contrast is very high (15 magnitudes for a 
600 K companion). The slope of the curves is different for each 
angular separation: it is steep for a separation of 1.5", and de- 
creases for increasing separations. The curves for separations of 
1.5", 1.0", and 0.5" reach the 80% level at different effective 
temperatures of 700 K, 900 K, and 1200 K, corresponding to 
contrasts of 13.6, 12.2, and 10.0 magnitudes respectively. This 
is consistent with the contrast values at which the 80% level is 
reached for the M0 star. 

Figure[9]summarizes the results of the quality factor for both 
M0 and GO stars at 10 pc at both LRS and MRS. It shows the 
effective temperature at which a value of q — 80% is reached as 
a function of angular separation. We clearly see the difference 
in regime between the two stars. A square law was fitted to the 
data to provide a general idea of the overall evolution, but more 
points are necessary to have a clearer idea of the real dependence 
on angular separation. 

For the M0 stars, the 80% level is reached over a small range 
in T e ff (less than 200 K). The influence of angular separation is 
then quite small, especially at MRS for which the curve is al- 
most flat. This is because at MRS we are limited by the sky level 
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Fig. 9. Effective temperature at which a value of 80% for the 
quality factor is reached in the case of MO and GO stars at 10 pc, 
as a function of angular separation, in LRS and MRS. A square 
law has been fitted on the points (dotted lines for LRS, dashed 
line for MRS) to show the overall evolution with angular separa- 
tion. The MRS results are only representative of the J band, and 
may change slightly for a larger wavelength range. 



Fig. 10. Influence of the planet mask size on the quality factor 
for 5 companions of increasing T e /f orbiting around an MO star 
at 10 pc and an angular separation of 1.0". Since the diameter 
of the planet PSF is wavelength dependent, the mask size is cal- 
culated in A/D units. The vertical dashed line marks the 2.5A/D 
limit that was used as default mask size in the simulations pre- 
sented in this work. 



for an M0 star at 10 pc. For GO stars, the influence of angular 
separation is far more significant: the range of Tgy covered to 
reach q = 80% is more than 600 K in LRS. At MRS, the range 
is only ~250 K, but compared to the evolution in M0 stars, this 
range is rather large. It finally shows that good characterization 
of companions with T e jf of 600 K orbiting at angular separations 
of 1 .0" around an M0 star at 10 pc can be achieved. For a GO star 
at 10 pc, the characterization of companions with a T e /f of 900 K 
is possible at separations of 1.0". 

4.2. Influence of varying parameters 

We study the influence of different data analysis parameters and 
procedures on the quality factor. To compare the different ef- 
fects, we focus on the case of an M0 star at 10 pc at LRS, with 
companions of various T e jf at an angular separation of 1.0". The 
MRS case was not considered given the fact that only 2 sets of 
simulation data in J band were available, but the conclusions for 
LRS would be applicable to MRS, with a possible scaling fac- 
tor related to the resolution difference. The two parameters that 



we studied were the planet mask size (Sect. 4.2.1 1 and the wave- 
length calibration (Sect. |4~2~2} . 

4.2.1. Planet mask size 



As explained in Sect. 2.2.2 a mask needs to be applied to the 
planet signal during the amplitude fitting of the model spectrum 
to the data to avoid any contribution of the companion signal to 
the level of scattered light. We studied the effect of changing the 
mask size as a function of the companion T e yy . The results are 



plotted in Fig. 10 



Two distinct regimes can be clearly identified: (i) the case of 
very faint companions (T e yj < 800 K) that cannot be recovered 
completely during data analysis because of noise and intrinsic 
limitations of the method, and (ii) the case of fairly bright or 
bright companions that can be recovered almost completely by 



the data analysis. In the first case, using a planet mask is unnec- 
essary because it either has no clear effect on the quality fac- 
tor (COND 700K model) or degrades the quality of the com- 
panion signal extraction (COND 600K model). For the 700 K 
case we see a very small improvement over the 0.0 - 2.5 A/D 
range, followed by a slow decrease towards larger sizes. For the 
600 K case, a slight amelioration of the quality factor is then 
followed by a clear degradation of the extraction when mask 
size increases. In this regime of high contrast where the noise 
becomes important and may become dominant over the signal, 
masking part of the signal only reduces the SNR and obviously 
degrades the data analysis. 

For relatively bright companions (T e jf > 800 K), we enter a 
second regime where the companion signal dominates the noise, 
and an almost complete recovery of the spectrum is achieveable. 
In this case, masking the planet signal is absolutely necessary. 
With a planet mask, the quality factors increase drastically for 
bright companions: with a mask of size 2.5A/D, a 6% gain is ob- 
served for the SETTL 1700K model. For colder companions, the 
effect of a mask of that size is less important, but not negligible. 
However, above 2.5A/D, the effect is very small. 

The origin of that limit is the encircled energy in the PSF 
of the companion: with a Strehl ratio of 90%, more than 77% 
of the energy is contained in a circle of diameter 2.5A/D. In our 
implementation of the method, we set the size of the planet mask 
to be this size. In the regime of faint companions, a degradation 
of ~1% or less for companions with T e ff < 800 K is accept- 
able compared to the gain of several percent that is obtained for 
brighter companions. In practice, it could be possible to decide 
whether or not a mask needs to be used on an individual case 
basis, but deciding which cases need a mask and which do not is 
not obvious a priori. This is why we decided to use a mask in all 
cases, without any variation as a function of T e ff . 



Arthur Vigan et al.: Exoplanet characterization with long slit spectroscopy 



9 



Systematic error @ 2.2 (nm) 




COND 600K 

COND 700K 

5ETTL 900K 

SETTL 1200K 

SETTL 1700K 



COND 600K 

COND 700K 

SETTL 900K 

SETTL 1200K 

SETTL 1 700K 



Systematic aberration error (nm) 
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Fig. 11. Influence of a systematic calibration error on the quality Fig. 12. Influence of an erroneous calibration slope on the qual- 
factor for 5 companions with increasing T e yj around an MO star ity factor for 5 companions with increasing T e yy around an MO 



l <ff 

at 10 pc and an angular separation of 1.0". For clarity, the error 
bars have not been plotted. For each model, they are almost con- 
stant: 4.0%, 0.7%, 0.25%, 0.2% and 0.2% for the COND 600K, 
COND 700K, SETTL 900K, SETTL 1200K, and SETTL 1700K 
models, respectively. 



4.2.2. Wavelength calibration 

Until now we assumed that the wavelength calibration is perfect 
i.e. that for each pixel column the corresponding wavelength is 
known precisely. Good wavelength calibration is a priori nec- 
essary to accurately rescale he spectrum and remove the wave- 
length dependence, and then be able to characterize the plan- 
etary companion. For the considered instrument, the calibra- 
tion of a long slit spectrograph with a resolution of R = 35 
can be achieved to a precision of ~0.2% from the slope of the 
wavelength dependence, and with a systematic error of less than 
15 nm. To understand the consequences of such errors in the cal- 
ibration more cleary, we simulated two different cases: (i) a case 
where a systematic error of a few nanometers was added to the 
calibration; and (ii) a case where the slope of the wavelength 
calibration was changed. 

Figure [TT] displays the effect of systematic errors in the cal- 
ibration for various companions around an M0 star at 10 pc and 
an angular separation of 1.0". Errors of 10 nm were introduced 
in steps of 0. 1 nm, and the resulting quality factor was plotted 
for various companion T e /f. The effect is globally a loss between 
2% and 3% in the quality factor, which is in any case almost neg- 
ligible. This does not appear to depend on the T e ff of the com- 
panion, although the effect is slightly worse for the COND 600K 
model. The slight loss is attributable to the correlation factor be- 
tween the input and output spectra that diminishes rapidly when 
the calibration error increases i.e. the different spectral features 
are no longer aligned, leading to errors in the characterization of 
the planet. For the COND 600K atmosphere model, the FWHM 
of the emission peak in the H band is ~108 nm, so a calibra- 
tion error of 10 nm produces a ~9% error in the position of the 
emission peak in the spectrum. 

For the second case, we introduced an error in the slope 
of the wavelength calibration. Figure 12] illustrates the result 
of introducing an error of 2% in the s ope, in steps of 0.05%. 
In the case that we are considering here, it means varying the 
slope from 13.328 nm/pixel to 13.872 nm/pixel in steps of 



eff 

star at 10 pc and an angular separation of 1.0". The bottom x- 
axis represents the simulated slope in nm/pixel and the top x-axis 
represents the resulting systematic error at 2.2 yt/m. 13.6 nm/pixel 
is the optimal slope. For clarity, the error bars have not been 
plotted. For each model, they are almost constant: 4.1%, 0.7%, 
0.25%, 0.23% and 0.27% for the COND 600K, COND 700K, 
SETTL 900K, SETTL 1200K, and SETTL 1700K models, re- 
spectively. 



68 x 10~ 4 nm/pixel. The results are similar to those for the sys- 
tematic error. The effect is globally the same for all the models, 
although the two coldest models (600 K and 700 K) appear to 
be slightly less sensitive to a systematic error than the warmer 
models. The COND 600K model undergoes a 1% loss for an 
error of 0.2 nm/pixel, while the SETTL 1700K model under- 
goes a 2% loss for the same error. However, for errors of 0.25% 
(34 x 10~ 3 nm/pixel) in the slope, the loss is smaller than 0.5% 
for all models. 

After analyzing the influence of two possible sources of er- 
ror in the wavelength calibration it appears that there is a clear 
necessity to achieve the smallest possible systematic error. An 
error of 10 nm (less than IAA) produces losses of several per- 
cent in the quality factor, even for bright companions. Reaching 
an uncertainty level smaller than the spectral interval AA is then 
an absolute requirement for a proper characterization of planets 
with our method. The influence of an error in the slope is also 
important, but the precision that can be achieved in the slope 
during the calibration procedure is small (0.25%). At this level 
of error, the loss in terms of quality factor is negligible. 



5. Conclusions 

We have implemented a promising method of characterizing 
planetary companions using long slit spectroscopy with coron- 
agraphy, which has been considered for very high contrast imag- 
ing instruments, such as VLT-SPHERE. The need to develop 
a specific method was mandatory to remove the scattered light 
residuals fully from the random speckles pattern in the spectra. 
Using the linear wavelength dependence of the speckle pattern, 
we were able to evaluate to high precision the precise contri- 
bution of the star to the spectrum, and remove this contribution 
before extracting a clean planetary companion spectrum. 
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The simulations, performed using IDL, allowed us to test our 
method on realistic data for various cases of contrast in the case 
of low (R = 35) and medium (R = 400) resolution spectroscopy 
with extreme AO and Lyot coronagraphy. The final gain of the 
method for a 1 hour exposure on MO and GO stars at 10 pc is 
of the order of 0.5 to 2 magnitudes in terms of contrast over 
the J, H and K bands, compared to the coronagraphic profile. 
Although not very high in K band, the gain is substantial in J and 
H, allowing us to study companions wi th T e ff as low as 600 K. 
Following the evolutionary models of Baraffe et al. (20031, a 
600 K companion in a 500 Myr system corresponds to a mass of 
~10M Jup . 

To estimate the method efficiency, a quality factor has been 
introduced to measure the correlation between the input plane- 
tary spectrum in the simulation and the extracted spectra, as well 
as the discrepancy between the two. This allowed us to be confi- 
dent that our method will allow precise characterization at LRS 
and MRS of companions as cool as 600 K around M0 stars at 
10 pc, or 900 K around GO stars at 10 pc, for angular separations 
larger than 1.0". For smaller angular separations, results should 
be extremely valuable for slightly warmer companions. 

Finally, we have estimated the influence of two key data 
analysis parameters at LRS. It is necessary to mask the planet 
signal during the analysis to avoid overestimating the scattered 
light level, which would lead to a mis-estimation of the plane- 
tary continuum. The method also relies on accurate wavelength 
calibration to be able to use the known wavelength dependence 
of the speckles to eliminate them. A systematic error as small as 
10 nm can have major consequences on the characterization of 
the observed objects. 
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